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SUMMARY 

This report details the work accomplished under NASA contract NAS2-6946, 
"Generation of NEP Heliocentric Trajectory Data. " The principal purposes of the 
study were to generate representative nuclear electric propulsion data for rendez- 
vous missions to the comet Encke using the variational calculus program HILTOP, 
to compare the data with equivalent data generated with the QUICKTOP program, 
and to propose approaches for storing and subsequently accessing the optimum tra- 
jectory and performance data in the QUICKLY program which is used for preliminary 
mission analysis studies. 

Optimum trajectories were generated with arrival at Encke 50 days prior 
to perihelion passage during the 1980 apparition. Four distinct trajectory classes 
covering a flight time range of 400-1400 days were investigated. The data pre- 

j, 2 

sented in this report consist of plots of the energy parameter J = J a dt, the pro- 
pulsion time ratio t /t_, and the arithmetic mean of the thrust acceleration 
i P f 

a = ( a Q a f) 2 as functions of the flight time for selected values of the launch excess 

speed in the range 0-8 km/sec. Tables of these three parameters and of the 

power to initial mass ratio p /m , the net spacecraft mass ratio m /m and 

o o no 

the jet exhaust speed c are also presented. The comparison of data indicates 
that the HILTOP and QUICKTOP program results for the mission investigated 
differ only minimally with the exception of the propulsion time. The final section 
of this report is devoted to a discussion of suggested modifications to the computer 
program QUICKLY to extend its capability and improve its accuracy. 
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INTRODUCTION 


In 1969, the Boeing Company, under contract to the Advanced Concepts 
and Missions Division (formerly Mission Analysis Division) at Moffett Field, 
California developed a low thrust trajectory optimization code CHEBYTOP^. 
Subsequently, ACMD personnel developed a driver package which calls CHEBY- 
TOP to obtain the optimum trajectory and then calculates vehicle masses and 
other parameters necessary for mission analysis studies. This program, known 
as QUICKTOP, was found to be easy to operate, relatively inexpensive to use, 
generally reliable in the sense of trajectory convergence and sufficiently accurate 
for most preliminary mission studies. In short, QUICKTOP has proved to be a 
useful and handy tool for low thrust mission analysis. Recently, more advanced 
versions of both CHEBYTOP^ and QUICKTOP^ were developed and are presently 
in widespread use. 

ACMD personnel also developed another low thrust mission analysis eom- 

(4) 

puter program known as QUICKLY . This program contains precomputed tra- 
jectory and performance data stored in polynomial form. Interpolation techniques 
are then used in conjunction with the characteristic length relations and other 
empirical formulas to estimate performance variations as a function of various 
propulsion system parameters. In generating data of selected missions for 
QUICKLY using QUICKTOP, considerable scatter, or erratic behavior, in the 
data was noted for certain parameters, principally the propulsion time, making 
the curve fitting of that data for QUICKLY particularly difficult. 

The study described in this report was undertaken for the following purposes 

(1) to generate optimal NEP heliocentric trajectory data using the variational 

(5) 

calculus program HILTOP v for a selected mission for which substantial scatter 
was observed; (2) to isolate any significant discrepancies between the results ob- 
tained from QUICKTOP and those from HILTOP; and (3) to consider alternate 
techniques in the program QUICKLY to expand its capability, improve its accuracy, 
and/or extend its flexibility while reducing the volume of data that must be stored. 
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NOTATION 


Parameter Definition 

a - semi -major axis; also, thrust acceleration 

a - arithmetic mean thrust acceleration 

a^ - final thrust acceleration 

a - initial thrust acceleration 
o 

b - coefficient employed in definition of efficiency factor 
c - jet exhaust speed of NEP system 
d - coefficient employed in definition of efficiency factor 
e - eccentricity 

f - thrust magnitude 

h - step function equal to zero if switch function is negative and equal to one 

if switch function is positive 

i - orbital inclination relative to ecliptic plane 

J - energy factor = Jj. ^ a 2 dt 

o 

- low thrust propellant tankage factor 
m - spacecraft mass 
m^ - net spacecraft mass 
m Q - initial spacecraft mass 
m ps " P r °P u ^ siLOn system mass 
p^ - beam power of NEP propulsion system 
P q - power into power conditioning system 

P q - optimum power level for a spacecraft using a specified launch vehicle 
opt 
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Parameter 


Definition 


p * - limiting power level below which more net mass can be achieved by off 
° loading the launch vehicle 

R - heliocentric position vector of spacecraft 

r - magnitude of R 

t - flight time 

t - initial time or launch date 

o 

t - accumulated time that propulsion system has operated during the mission 

v^ - hyperbolic excess speed 

a - specific propulsion system mass 

77 - total efficiency of the propulsion system 

77' - derivative of efficiency with respect to jet exhaust speed 

A - primer vector 

X - magnitude of A 

X o - multiplier adjoint to initial thrust acceleration 

X - multiplier adjoint to jet exhaust speed 
c 

X^ - multiplier adjoint to mass ratio 

/i - gravitational constant of the sun 

v - instantaneous mass to initial mass ratio 

v - net spacecraft mass ratio 
n 

V - low thrust propellant mass ratio 
P 

- NEP propulsion system mass ratio 
v - propellant tankage mass ratio 
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Parameter 


Definition 


CT - engine switching function 

r - perihelion passage date 

co - argument of perihelion 

Cl - ecliptic longitude of ascending node 
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DATA GENERATION 


1. Study Guidelines. 

The mission selected for study in this contract is the comet Encke 
rendezvous mission, arriving at the target 50 days prior to perihelion passage 
in the 1980 apparition. This corresponds to a Julian date of arrival of 2444530 
(October 17, 1980). The study was restricted to the heliocentric portion of the 
mission with the arrival excess speed fixed at zero. Launch excess speeds in 
the range 0-8 km/sec were investigated for total flight times in the range 
400-1400 days. A launch vehicle independent formulation in which all masses 
are referenced to the initial spacecraft mass was employed. The power level 
and jet exhaust speed were optimized to yield maximum net spacecraft to initial 
mass ratio. Data were generated for a baseline specific propulsion system mass 
of 5 kg/kw with sensitivity data cases being run for specific masses up to 40 
kg/kw. A propellant tankage factor of 0. 03 was assumed throughout. Compari- 
son data obtained from the QUICKTOP program were provided by the Technical 
Monitor. 

2. Formulation . 

The optimal trajectory and performance data for this study were ob- 
tained with the Heliocentric Interplanetary Low-Thrust Trajectory Optimization 
(5) 

Program (HILTOP) . This program utilizes the indirect variational calculus 
method for optimizing a three-dimensional low thrust heliocentric trajectory 
between assumed massless bodies moving in prescribed orbits about the sun. 

An analytic ephemeris is employed which yields orbital positions in a helio- 
centric ecliptic frame of date. The orbital elements of Earth are expressed as 
quadratic functions of Julian century in this frame. The assumed constant ele- 
ments for the comet Encke in the 1980 apparition are as follows: 

a = 2.218 AU e = 0.847 i = 11.95 deg 

60 = 185. 98 deg C= 334.19 deg T = 2444580 (Dec. 6, 1980) 
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A trajectory is obtained by numerical integration of the equations of 


motion 


R = h — R 

cr v X 3 
r 


v = - h 


a c 


simultaneously with the adjoint equations 


A= ^ (R- A)R--^ A 

5 o 

r r 


X =h — n X. 
y a 

The step function h^ is chosen on the basis of the sign of the switch function 

CT = X - — X as follows : 
c v 


r 0 if cr<0 
a 1 1 if o > 0 


The initial position used to start the integration is that of Earth on the specified 
launch date, while the initial velocity is set equal to the velocity of Earth plus 
the specified hyperbolic excess speed applied in the direction of the primer 
vector, A. The initial mass ratio v is, of course, equal to one, and the initial 
values of A , A , and X^ are all guessed. The integration is terminated upon 
reaching the desired flight time. 

An optimum trajectory for a specified trajectory class, flight time, and 

launch excess speed is obtained by solving the boundary value problem for a , c, 

. ^ 

and the initial values of A, A , and X^ to yield the desired final position and 
velocity, and to satisfy the three transversality conditions 


X 

X 


,<V = 1+k t 


(t ) = m /a 
a f ps o 
o 
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\ (t ) = - m ( i-i') 

C ' f PS \ C T) / 


where X and X are obtained by numerical integration of the equations 
a c 

o 

X = -h - 
a. o v 
o 

a 

X* = - h X . 
c a 2 v 
c 


The initial values of both X and X are zero. Of course, the desired final 

a c . 

o 

position and velocity is taken to be that of comet Encke on the specified arrival 
date. 


The powerplant is assumed to deliver a constant power p to the power 

o 

conditioner. A single efficiency factor, rj, is employed to relate beam power 
to the input power as follows : 


p b “ " P o • 


with 77 assumed to be a function of the jet exhaust speed. 

2 




b c 

2 2 
c +d 


The values of the constant coefficients b and d used in this study were 0.842 and 
16 km/sec, respectively. The thrusters are assumed to operate at constant jet 
exhaust speed and mass flow rate throughout all propulsive maneuvers. The 
dynamic relationships that apply during thruster operations are: 


1 . 2 

p b = “ 2 m 0 =7,p o 


f = - me = a m . 

o o 
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Since the mass flow rate is assumed constant the spacecraft mass and mass 
ratio are linearly decreasing functions of propulsion time. 


The spacecraft is assumed to consist of propulsion system, propellant 
and tankage. Any additional mass available is termed net spacecraft mass. 
Thus, in terms of mass ratios 


v +v+v+v=l 

ps p t n 


where 


v =ap /m =aa c/2 ri 
ps o o o 


V ’ = k .v 

t t p 


Hence, 

v = 1 - aa c/2T7 - (1 +kj v , 
n o t p 

which is to be maximized in the optimization procedure. 

Other relationships that are necessary for comparison of results with 
those of the QUICKTOP program are the energy factor 


J = f 4 f a 2 dt = a c v /(1-v ) , 

W n TA TY 


o P 


and the geometric mean thrust acceleration 


a = (a a ) 2 = a / J\-v~ . 
of o p 


A useful relationship that results from the constant mass flow rate yields the 
propulsion time as a function of J and a: 

t - J/a 2 . 

P 
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3. 


Results. 


Optimal trajectories for the Encke rendezvous mission in the flight time 
range of 400-1400 days may be divided into four distinct trajectory classes. 
These different classes arise because the same geometric configuration for 
transfers from Earth to the position Encke occupies on the desired arrival 
date repeats annually. Given a trajectory with a specified flight time, there 
will exist another trajectory with exactly the same travel angle and arrival date 
but with a flight time exactly one year longer. The longer flight time of one 
year is achieved by simply increasing the aphelion distance through which the 
trajectory passes. The sketch in Figure 1 depicts trajectories in each of the 
four classes that were considered in this study. The sketch depicts four tra- 
jectories arriving at Encke on the same arrival date and departing Earth at the 
same heliocentric location, but on dates exactly one year apart. The earlier 
launch dates (and hence the longer flight time) correspond to the larger aphelion 
distance. Hereafter, the terms short, intermediate, long and extra long flight 
time class will be used to differentiate among the four trajectory classes. 


For each of the flight time classes, tables are presented which list the 

optimal values of six parameters as a function of launch excess speed and flight 

time. The six parameters include the energy parameter J, the propulsion time 

to flight time ratio t /t^, the mean thrust acceleration a, the power to mass ratio 

p /m , the net spacecraft mass ratio m /m (=u ), and the jet exhaust speed c. 
o o non 

The values presented are those obtained from the HILTOP program. When avail- 
able a corresponding data point obtained from the QUICKTOP program is included 
in parentheses. 


In Table 1 are presented the data for the short flight time class. The 

table is divided into six parts which are labeled la through If and which contain 

the parameters J, t /t_, a, p /m , V , and c, respectively. Each parameter 
p t o o n 

is tabulated as a function of launch excess speeds at increments of 1 km/sec over 



the range 0-8 km/sec and as a function of flight time at 50 day increments over 
the range 400-600 days. The only comparison points available for this class of 
solutions were at a flight time of 400 days. This comparison indicates relatively 
minor discrepancies between the two methods of solution, the largest being in the 
energy factor J which varies from 5-10 percent. It is particularly interesting to 
note that this discrepancy in J is not matched by a corresponding discrepancy in 
the net spacecraft mass ratio. This is easily explained as follows. Letting 6 J 
denote a discrepancy in J and 61 / a discrepancy in the propellant mass ratio V 
then one would expect that 




P 


E 

aj 


6 J 


and since 


v 

P 


J 

a c + J * 
o 


then, neglecting any discrepancies in a^ and c, and noting that 
obtains 


6 v — - 6 v one 
n p 


5v = - 
n 


a c 
o 


(a c + J) 
' o 


6 J . 


2 3 

For the problem at handle and J are about equal and of magnitude 50 m /sec . 
Consequently, 61 / should be about two orders of magnitude less than 6 j, which 
is close to the difference exhibited in the tables. 


In Table 2 are presented the data for the intermediate flight time class. 
The organization of the data is the same as for Table 1 with flight times at 50 
day increments in the range 550-700 days and at 100 day increments thereafter 
to 900 days. The grid points for this class which yield direct comparisons of 
the HILTOP and QUICKTOP results indicate differences of 5-10 percent in J, a, 
and c and of about 1-4 percent in net spacecraft mass ratio. A certain con- 
sistency is evident in the comparisons in that the QUICKTOP results are always 
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lower in J and a and higher in v and c than the corresponding results of HILTOP. 
The most obvious discrepancy in the results of the two programs is in the pro- 
pulsion time ratio. Whereas the QUICKTOP program indicates continuous pro- 
pulsion to be optimum in all but one of the 14 grid points shown, the HILTOP 
program yields propulsion time ratios varying from about 0.6 to 0.9. Clearly, 
since the QUICKTOP yields more favorable performance estimates on the basis 
of optimal propulsion time, the difference in the output of the two programs would 
widen if HILTOP were constrained to continuous thrust solutions. For example, 
a few typical cases run with HILTOP indicated that forced continuous thrust re- 
sulted in a 2-5 percent increase in J and 0.5-1 percent decrease in v as 

n 

compared to optimal propulsion time results. The cases considered included 
flight times of 600, 700, and 800 days at a launch excess speed of zero. The 
optimal propulsion time ratios for these cases were about 0.85 (see Table 2b). 

The long flight time class data are presented in Table 3 at 100 day in- 
crements in flight time over the range 900-1300 days. The results of the two 
programs compare more closely for this trajectory class than for the intermediate 
flight time class. With few exceptions, the results for this class are within two 
percent. Differences in propulsion time ratio, if they exist, are not so noticeable 
in this case because both programs yield continuous propulsion solutions at most 
of the grid points available for comparison. Similarly in Table 4, the single 
point available from QUICKTOP for the extra long flight time class compares 
very favorably with that from HILTOP. A difference of about two percent may be 
noted in J and c. The vacant grid points in Tables 1-4 are due to numerical pro- 
blems experienced with HILTOP when attempting to attain those specific optimum 
trajectories. 

To facilitate an understanding of the behavior and trends of the data pre- 
sented in the tables, plots were prepared of the parameters J, a and t /t. as 

P 1 

functions of both the flight time and the launch excess speed. The three para- 
meters are shown in Figures 2-4, respectively, as functions of flight time for 
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curves of selected launch excess speeds. These figures were prepared with 
the flight time plotted on a linear rather than the logarithmic scale suggested 
in the Statement of Work. The reasons for this choice were that launch 
opportunities repeat annually, which suggests a linear scale, and that compressing 
data at the longer flight times appeared to amplify rather than reduce problems 
one would have in curve fitting the data. 

The results in Figures 2-4 are not particularly favorable from the stand- 
point of curve fitting for the QUICKLY program. One will note for both J and a 
that the general shapes of the curves change drastically from one trajectory class 
to the next and, to a lesser extent, from one excess speed to another. This 
would suggest potential difficulty in finding an acceptable general form of an equa- 
tion to which the data may be accurately curve fitted. Figure 4 showing the propulsion 
time ratio has even more scatter. Here, with one exception, no attempt was made 
to fare curves through the data points. The linear dashed lines between data points 
are included to facilitate distinguishing the class of solutions to which each point 
belongs. Sufficient detail was developed to define the shape of the curve of pro- 
pulsion time ratio versus flight time for an excess speed of 8 km/sec in the long 
flight time class and the results are shown by the solid curve in the Figure 4. 

The curve exhibits six distinct arcs within the flight time range 900-1200 days. 

From 900-920 days, the ratio changes very rapidly from about 0. 73 to 1. 0; from 
920-940 the solution indicates continuous propulsion; from 940-980 the ratio is 
parabolic in shape, exhibiting a minimum value of about 0.85; and from 980-1200 
there are three nearly linear segments with abrupt changes in slope at about 1080 
and 1125 days. This rather unusual and unexpected shape was deduced from ob- 
servations of the switch function time histories associated with the eight available 
data points along the curve. These time histories are shown in Figure 5. 

To understand Figure 5, one must first realize that the switch function O' 
is used to optimally switch the engine on or off when it passes through zero. When 
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a < 0 , the engine is off and, conversely, when cr > 0, the engine is on. This 
figure is useful in understanding the propulsion time ratio curve in Figure 4 in 
that it shows the duration and location of individual coast phases in each solution. 
For example, the 900 day solution exhibits a negative switch function for about 
240 days starting at Julian date 4130. At all other points, the switch function is 
positive. The minimum on the 910 day curve is shifted upward and to the right 
such that the coast phase is reduced to about 180 days. Otherwise, the curve 
closely resembles that for the 900 day case. Proceeding to the 920 day solution, 
one will note that the minimum has shifted further to the right and also upward 
to the extent that the function is always positive which yields the propulsion time 
ratio of one as seen in Figure 4. Meanwhile, with increasing flight time, a mini- 
mum of the switch function has developed early in the flights, around the Julian 
date 3700 for the 920 day case. For the 950 day case, it is seen that the late 
minimum has all but disappeared while the early minimum is noticeably more 
pronounced and has dropped below zero. Thus, the continuous thrust solutions 
are only optimal for a small flight time range of about 920-940 days. As the 
flight time increases above 940 days, the early minimum dips well below zero, 
resulting in a sizeable coast phase, and then quickly returns to near zero to give 
the parabolic dip shown in Figure 4. The minimum then hugs the zero line as is 
seen in the Figure 5 for flight times of 1000, 1050 and 1100 days. The minimum 
raises above the zero axis for the 1200 day solution, thereby eliminating that 
coast phase. This is the source of the abrupt change in slope of the curve in 
Figure 4 at a flight time of about 1125 days. The change in slope at about 1080 
days is due to the new coast phase that is introduced in the solutions at launch 
and is first evident in the curves in Figure 5 for the 1100 day case. Thus, the 
unusual shape of the propulsion time curves as a function of flight time are seen 
to result from the rapid changes in the switch function time histories as the flight 
time is varied. If the curves of Figure 5 were drawn in perspective with flight 
time as the third coordinate, one would observe a rather complex surface re- 
sembling a rough terrain with high mountains and deep valleys. With this picture 
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in mind it then becomes somewhat easier to accept the erratic behavior of the 
propulsion time curves shown in Figure 4. The lesson to be learned from this 
example is that the propulsion time ratio must be considered (at least for the 
Encke rendezvous mission) as a rapidly fluctuating variable that does not lend 
itself well to a curve fitting approach. 


The parameters J, a and t /t^ are plotted as functions of launch excess 

speed in Figures 6-8, respectively. These curves are seen to be quite smooth 

functions of the independent parameter and should be somewhat more amenable 

to curve fitting than the corresponding curves plotted as a function of flight time. 

The parameters J and a for long flight times within each class appear to match 

quite favorably a linear curve fit on the semi-log grid over the range of launch 

v ’s shown. The shorter flight times in each class (e.g. , 900 days-long and 
00 

1300 days-extra long) will probably require a higher order fit. It is particularly 
interesting to note the striking similarity in the shapes of the curves for J and a. 
It would appear that identical functional forms could be used for the two variables. 
Because of this it is probably preferable to avoid curve fitting the propulsion time 
ratio since it may be immediately obtained from J and a. 


The sensitivities of J, a and t /t, to the specific propulsion system mass 

P f 

a were obtained for several selected cases in the intermediate flight time class. 
For a 700 day mission solutions were obtained for several values of a in the 
range 5-40 kg/kw and for launch excess speeds from 0-8 km/sec. In addition, 
cases were obtained for the same interval in a holding the launch excess speed 
fixed at 4 km/sec and varying the flight time in the range 600-800 days. The re- 
sulting variations in the parameters J, a and t /t are presented graphically in 

P ^ 

Figures 9-11, respectively. Although somewhat non-linear, these curves are 
very smooth functions of a, implying little difficulty in curve fitting the data if 
desired. The smoothness in the data also indicates that choosing a different 
nominal value of a for the generation of the basic data in QUICKLY would not 
influence the degree of difficulty in curve fitting data as a function of flight time. 
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4. 


Concluding Remarks. 

From the data generated for the Encke rendezvous mission in this study, 


it is possible to draw the following conclusions regarding the curve fitting of the 
data in QUICKLY: 


1. Curve-fits as a function of flight time must be restricted to re- 
latively short arcs within each class of solutions if accurate propulsion time 
estimates are required. Concomitant with this requirement will be the need 
for initially generating the optimal trajectory data on a small flight time grid 
(~ 10 days). Performance estimates without regard to propulsion time may 
be obtained on a somewhat larger grid (~ 50 days). Depending upon the accuracy 
desired, it maybe necessary to choose different functional forms for the curve 
fit for the different classes of solutions and possibly for the different excess 
speeds within a class. 


2. Of the three parameters J, a and t^/t^, only two are independent, 
and since J and a appear to be smoother functions of both flight time and launch 
excess speed, it is recommended that they be chosen for curve fitting. 


3. Curve fits of J and a as functions of launch excess speed or specific 
propulsion system mass for a given flight time may be adequately achieved with 
low order equations that apply over the entire ranges of v o and a investigated in 
this study. In many cases a linear fit on a semi-log grid is appropriate. 


A comparison of the data generated in this study with that generated with 
the QUICKTOP program and supplied by the Technical Monitor leads to the follow- 
ing conclusions for the specific mission investigated: 

1. The net spacecraft mass ratios obtained from the QUICKTOP and 
HILTOP programs for the same case compare very favorably, usually differing 
by less than one percent. 
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2. Parameters such as J, a and c usually compare within 10 per- 
cent and frequently within 5 percent. 

3. Propulsion time ratio appeared to be the parameter in which the 
greatest discrepancies existed in the results of the two programs. The largest 
discrepancies occurred for the intermediate flight time class for which QUICKTOP 
generally predicted continuous thrust whereas HILTOP predicted propulsion times 
of 60-90 percent of the flight time. 
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CONSIDERATIONS OF QUICKLY 


A second, and essentially independent, task of this study was to review 
the design of the program QUICKLY and to suggest changes to the program which 
would extend its capability and/or improve its accuracy while retaining its de- 
sirable attributes of simplicity and speed of operation. Briefly, the manner in 
which QUICKLY currently operates, as described in Reference 4, is as follows. 
Stored in data arrays are coefficients of specific polynomials which yield the 

performance parameters J and t as a function of the flight time and the sum 

P 

of the launch and arrival excess speeds (launch excess speed only for flyby 
missions) for specific missions of interest. These coefficients were determined 
independently from least squares curve fits of precomputed optimal trajectory 
data generated with open power and jet exhaust speed to yield maximum net space- 
craft mass ratio for specified values of a, and any other structural factors. 

Given J and t for specified flight time and set of hyperbolic excess speeds, 

P 

analytic equations are then entered which allow one to compute the optimum values 
of the spacecraft parameters such as c, a Q , V , a » etc., in closed form. 

Thus, it is not necessary to curve fit these spacecraft parameters as they are 
not independent. Empirical relations are then employed to obtain estimates of the 
changes in performance and spacecraft parameters due to changes in design para- 
meters, such as a, and due to the imposition of constraints, such as fixed power 
and/or limited propulsion time. The empirical equations used to handle the con- 

(CK \ 

straints involve variations of Zola's characteristic length' ’ method which assumes 
invariance of the characteristic length for a specific mission over a wide range of 
system parameters. The characteristic length is, roughly, the second integral of 
the thrust acceleration. 

Except for the extent to which precomputed mission data is stored and 
available in the program, the greatest limitations in QUICKLY would appear to 
come from two sources: (1) the accuracy to which the precomputed data are 
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curve fitted, and (2) the validity of the empirical relations employed. Each of 
these sources will now be considered with specific attention being given to the 
numerical results obtained for the Encke mission described in the preceding 
section. One additional factor which has nothing to do with the design of 
QUICKLY but which would appear to represent a limitation in QUICKLY is the 
accuracy of the precomputed data itself. 

The range of flight times and number of classes of solutions considered 
in this study for the Encke rendezvous mission undoubtedly exceed those for 
which QUICKLY was originally envisioned to handle. This supposition combined 
with the fact that the Encke mission requires relatively large changes in inclina- 
tion and eccentricity should be sufficient reason to suspect that the optimal per- 
formance data for the Encke mission may not be as well behaved as, say, a 
mission to one of the planets, and hence may present more severe problems in 
curve fitting the data. This is indeed the case as is obvious in Figures 2-4. 

Little can be stated conclusively in this report concerning explicit curve fitting 
techniques to adequately handle the Encke rendezvous mission because the scope 
of work did not permit sufficient in-depth study to cover this aspect of the pro- 
blem. However, due to the striking similarity of the curves of J and a, as is 
seen in Figures 2, 3, 6 and 7, and due to the erratic behavior of the propulsion 
time ratio curves, a study of techniques for curve fitting J and a rather than 

J and t seems warranted. Such a study should include an examination of the 
P 

resultant error in t and the implications of the error as it affects the charac- 

P 

teristic length relations and appropriate procedures to follow when the calculated 

t is greater than t . 

P f 

The use of empirical equations to describe relations among a set of variables 
carries with it a certain amount of risk in that one can only be certain that it applies 
over the ranges of the variables that were used in deriving the equations. This is 
particularly true if the functional form of the relationship was derived from a num- 
erical rather than an analytic basis . If the assumed form of the equation incorrectly 
describes the relationship, erroneous results will obviously be obtained from the use 
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of the empirical equation. A pertinent example of this is found in the data presented 
in the preceding section. In Figure 9, the behavior of J as a function of a is 
shown for five different flight time/launch excess speed combinations. In four of 
the five cases the slopes of the curves at a = 5 are negative. The empirical 
equation appearing in QUICKLY which is intended to describe the variation in J 
as 0! is varied will not account for a change of sign of slope. 

The use of the characteristic length equations also represents an empirical 
approach. Furthermore, a number of additional empirical relations were developed 
for QUICKLY, again based on numerical observations, which attempt to correct the 
characteristic length to account for various mission constraints. Because of a lack 
of knowledge as to the scope of the cases observed in deriving these corrections, 
there is some question regarding the general applicability of the equations. A 
method is available that appears to eliminate the need for empirical scaling to 
cover constrained power cases or variations in a. The method also has other 
attributes as well as some disadvantages. The remainder of this report is devoted 
to a discussion of this method. 

(7) 

Recently a concept was published which describes an approach and a 
set of assumptions under which one may generate optimal trajectory data that is 
independent of the launch vehicle, the power level, the specific powerplant mass 
and the propulsion system efficiency. Here the term independent implies that the 
optimum trajectory may be defined without regard to any of the four parameters 
and therefore effects of variations of any or all of the parameters on the perfor- 
mance may be determined exactly without regard to any approximate or empirical 
formulas. These features are achievable by redefining what is meant by an 
optimal trajectory. The optimal data discussed in the preceding section was gen- 
erated under the criterion of maximizing the net mass ratio with open power and 
jet exhaust speed and specified launch and arrival dates and hyperbolic excess 
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speeds. The independence cited above is achieved by declaring the performance 

index to be the net mass to power ratio m /p rather than the net mass ratio. 

no 

Although this may seem on the surface to be a somewhat contrived requirement, 
it actually proves to be the best choice of the performance index for some 
cases of interest - specifically those for which the power level is constrained 
to a level well below that which yields the true maximum net mass for a given 
launch vehicle. This is seen to be true from the following discussion. 

Suppose a mission to a given target is specified in terms of the launch and 

arrival dates and hyperbolic excess speeds, and several trajectories for a range 

of values of the power p are optimized to yield maximum net spacecraft mass. 

o 

We will assume that the same initial mass, corresponding to the payload capabi- 
lity of a selected launch vehicle at the specified launch excess speed, is used at 
all power levels, and also assume that the jet exhaust speed is optimized at each 
power level. Then a typical plot of the maximum net spacecraft mass as a function 
of power level will appear as shown by the solid curve in the sketch of Figure 12. 

The absolute maximum of m occurs at a power level of p whereas the maxi- 

n °opt 

mum value of the ratio m /p occurs at a power level p * corresponding to the 

no o 

tangency point of the dashed line extending from the origin and the solid curve. An 
interesting property of optimal solutions with constrained power was noted inde- 
pendently by Meissinger^ and Zola^. This property may be observed in Figure 
12 by noting that at power levels below p^*, a larger net spacecraft mass can be 
achieved by following the linear dashed line than by following the solid curve. To 
follow the dashed line requires that one employ the jet exhaust speed of the tangency 
point and scale proportionately the power level and all mass components. If this 
is done, the trajectory followed will be the same for any point on the dashed curve. 
The reduced initial mass at any point on the dashed curve represents the optimum 
initial mass for the corresponding power level. Consequently, if a larger initial 
mass is used, the additional propellant and tankage requirements necessary to 
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meet the specified end conditions will exceed the increase in initial mass thereby 

causing a net reduction in the net spacecraft mass. Now, note that the tangency 

point of the dashed and solid curves represents the maximum of the ratio m /p 

n o 

that is achieved anywhere on the solid curve. It then follows that the choice of 

this ratio as the performance index will indeed yield the solution for maximum 

net spacecraft mass achievable for powers less than p *. 

o 

To better understand how the approach is implemented in a mission study, 

suppose a set of optimal trajectories were obtained for maximum m /p over a 

n o 

range of launch excess speeds. For each value of launch excess speed there would 
result from the optimization corresponding optimal values of a Q and c as well 
as the propellant fraction v and the propulsion time t^. Of course, the para- 
meters J and a are also available. Note that each of these parameters is invariant 
along the dashed curve of Figure 12. From the mass relations given in the preced- 
ing section for the Encke rendezvous mission one may write 


m 


n 


a c 
o 


[l - (l +k ) v ]-a 
t p 


m 

o 27? 


p a c 
o o 

Now, due to the dependence of a , c and v on launch v , both m /p and 

o p “no 

m /p are also functions of v . Suppose one now selects values of a, 77 and 
o o “ 

p and plots m as a function of launch v and selects the case (i.e. , the v ) 
o n “ oo 7 

at which m is maximized. But from the equation for m /p it is clear that 
n no 

although the choice of a affects the value of m^, it does not affect the value of 
v^ at which m^ is maximized. Similarly, a change in the value of p^ or 7? 
will affect only the magnitude of m n and not the location of the maximum. Con- 
sequently one accounts exactly for any changes in these parameters through the 
above equations. Furthermore, to consider a specific launch vehicle, one simply 
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multiplies the function m Q /p o (=2 17/a c) by the prescribed power. This, in 

effect, yields a curve of optimum initial mass as a function of launch v . The 

intersection of this curve with that of the payload capability of the desired launch 

vehicle, plotted on the same scale, yields the appropriate launch excess speed 

and, in turn, the maximum net mass capability for the specified launch vehicle 

and power level. A maximum power level that one should attempt to use with 

this approach for a given launch vehicle may be determined by dividing the pay- 

load capability of the launch vehicle at a given launch v ro to obtain a curve of the 

optimum power level p * for that launch vehicle as a function of launch v . 

0 00 

The peak of this curve represents a limit to the power level beyond which the 
optimum initial mass always exceeds the payload capability of the launch vehicle 
and hence no solution exists. 

To further clarify the use of this multi -parameter independent mode 
(MPIM) a set of optimal data for a specific mission were generated. The mission 
chosen was the 700 day Encke rendezvous in the intermediate flight time class 
for which data are also available in the Table 2 for maximum net mass ratio. The 
performance data generated using the new mode are tabulated as a function of 
launch excess speed and are presented in Table 5. The parameters include all 
those presented in Table 2 for the same mission so that a direct comparison of 
the results of the two modes can be made. Such a comparison leads to the follow- 
ing observations: 

(1) The energy factor J differs only slightly in the two approaches. 

(2) The propulsion is continuous for the MPIM. This will not always 
be the case. However, longer propulsion times are expected to result in general 
using the MPIM. 

(3) The power to mass ratio p^/m^ is considerably smaller (by 

about a factor of 5) for the MPIM. Of course, in Table 5 p /m represents 

o o 

p */m whereas in Table 2 it represents p /m . 

00 o^o 

opt 
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(4) 


Similarly, because of the different power levels at which the 


points apply, the net mass ratio is much smaller for the MPIM, 

(5) Possibly the most striking difference in the results of the two 
methods is in the jet exhaust speed c. Whereas the optimum value using the 
original method was over 120 km/sec, it is reduced to about 26 km/sec using 


the MPIM. This is because the reduction in power from p 


to p * is achieved 
o , o 
opt 


primarily by reducing the jet exhaust speed and only secondarily by reducing the 
thrust acceleration. 


Data for the 700 day mission in Table 2 were employed to develop the NEP 
performance data for the Titan III D/Centaur launch vehicle. The data are pre- 
sented in Figure 13 and include the initial NEP spacecraft mass (assumed equal 

to the payload capability of the launch vehicle), the optimum power p (equal 

opt 

to the ratio, p /m from Table 2 multiplied by the initial mass), and the net space- 
o o 

craft mass (equal to the product of v ^ and mj. The parameters are plotted as a 
function of launch excess speed. The peak of the net mass curve corresponds to 
the solution generally obtained for this launch vehicle under the assumption of 
optimum power and jet exhaust speed. Note that if constrained power data is de- 
sired, it is essential to use the characteristic length technique to scale this data. 
Similarly empirical equations are required to approximate variations in these para- 
meters due to changes in a or T). The data from Table 5 were used in developing 
corresponding curves in Figure 14 for the same launch vehicle implementing the 
MPIM. Plotted as a function of launch excess speed are the launch vehicle pay- 
load capability, the power level P q * for this launch vehicle, and the net mass 
associated with the power level P q *. In addition, dashed curves are presented 
for the case of power fixed at 15 kilowatts. Shown are the optimum initial mass 
for this power level and the net spacecraft mass associated with that solution. 

Note that these dashed lines intersect the corresponding solid lines and terminate 
at the launch excess speed for which p * = 15 kilowatts. For all lesser excess 
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speeds the optimum initial mass is less than the payload capability of the launch 
vehicle, and for those cases the launch vehicle could be off-loaded. To consider 
a different a in this example, all that is required is to evaluate the difference in 
propulsion system weight and transfer to or from the net mass on a one-to-one 
basis. Similarly, a change in efficiency is accommodated by scaling the initial mass 
and the sum of the propulsion system and net masses by the ratio of the new to the 
old efficiency factor. Then subtracting the propulsion system mass from the 
scaled sum yields the new spacecraft mass. 

The curve fitting of parameters for the MPIM should prove no more diffi- 
cult than the original approach. No change in the choice of variables curve fitted 
is apparent at this time although that possibility should be considered. It will be 
necessary to re-derive the appropriate expression for optimal c, given J and t 
or a, since the performance index has changed. An alternative is to generate 
optimum trajectory data for specified rather than optimal values of c and either 
curve fit the data as a function of c or use the characteristic length relations to 
approximate effects of changes in c. Similar alternatives exist with respect to 
constrained propulsion time. Finally, it is important to note that the MPIM is 
adaptable to other types of missions, such as orbiters and flybys, with no more 
difficulty than with the original approach. The principal limitation of the MPIM 
is that power levels should not exceed p *. 
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Table la - Energy Factor, J (m /sec ) 
Short Flight Time Class 


Launch v 

03 


Flight 

Time (days) 



400 

450 

500 

550 

600 

0 

54.55 (50.3) 

46.05 

45.58 

44.96 

41.97 

l 

50.76 

44.51 

44.84 

43.71 • 

40.12 

2 

47.40 

43.36 

44.37 

42.71 

38.60 

3 

44.42 

42.44 

44.11 

41.97 

37.31 

4 

41.82 (41.4) 

41.63 


41.42 

36.19 

5 

39.56 

40.97 


41.00 

35.20 

6 

37.61 (41.4) 

40.45 


40.70 

34.33 

7 

35.94 

40.08 


40.51 

33.56 

8 

34.54 (31.0) 

39.85 


40.41 

32.90 


Table lb - Propulsion Time Ratio, t /t 
Short Flight Time Class 1 


Launch v 


Flight Time (days) 








0 

.7596 (.73) 

.7633 

.7180 

.6915 

.6946 

1 

.7610 

.7597 

.7078 

.6829 

.6908 

2 

.7599 

.7572 

.6997 

.6771 

.6870 

3 

.7565 

.7373 

.6942 

.6887 

.6830 

4 

.7515 (.72) 

.7182 


.6991 

.6788 

5 

.7458 

.7003 


.7089 

.6745 

6 

.7400 (.72) 

.6840 


.7184 

.6700 

7 

.7347 

.6695 


.7273 

.6722 

8 

.7304 (.72) 

.6566 


.7357 

.6757 
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Table lc - Mean Thrust Acceleration, a, (10 m/sec ) 
Short Flight Time Class 


Launch v 

00 


Flight 

Time (days' 

) 


400 

450 

500 

550 

600 

0 

14.41 

12.46 

wmm 

W 

10.80 

1 

13.89 

12.28 



10.58 

2 

13.43 

12.14 

12.12 


10.41 

3 

13.04 

12.17 

12.13 


10.27 

4 

12.69 

12.21 


11.17 

10.14 

5 

12.39 

12.27 


11.03 

10.03 

6 

12.13 

12.33 


10.92 

9.94 

7 

11.90 

12.41 


10.83 

9.81 

8 

11.70 

12.49 


10.75 

9.69 































Table Id - Power to Mass Ratio, p /m (kw/kg) 
Short Flight Time Class 


Flight Time (days) 




.0518 

.0490 


.0481 

■ 

.0508 

.0484 

KSifl 

.0480- 


.0498 

.0479 

.0483 

.0478 


.0489 

.0477 

.0483 

.0476 

.0461 

.0480 

.0476 


.0475 

.0457 

.0472 

.0475 


.0473 

.0454 

.0466 

.0474 


.0473 

.0451 

.0460 

.0474 


.0472 

.0448 

.0454 

.0474 


.0472 

.0446 



Table le - Net Spacecraft Mass Ratio, m /m 
Short Flight Time Class n ° 


Flight Time (days) 


400 

450 

500 

550 

600 

3337 (.335) 

.3768 

.3795 

.3834 

.4006 

.3517 

.3849 

.3834 

.3901 

.4112 

3685 

.3911 

.3858 

.3956 

.4202 

,3842 

.3960 

.3872 

.3999 

.4280 

3986 (.399) 

.4003 


.4032 

.4350 

,4117 

.4039 


.4056 

.4412 

4234 (.425) 

.4067 


.4075 

.4469 

,4338 

.4087 


.4087 

.4520 

4428 (.443) 

.4099 


.4094 

.4565 


Table If - Optimum Jet Exhaust Speed, c (km/sec) 
Short Flight Time Class 



Flight Time (days 


400 

450 

500 

550 

600 

74.45 (73) 

80.08 

81.25 

83.94 

88.92 

74.95 

79.90 

81.09 

83.89 

89.19 

75.32 

79.73 

80.88 

83.95 

89.41 

75.55 

78.95 

80.72 

84.92 

89.57 

75.70 (74) 

78.18 


85.74 

89.67 

75.78 

77.43 


86.48 

89.72 

75.84 (75) 

76.72 


87.17 

89.73 

75.90 

76.05 


87.81 

90.10 

75.99 (76) 

75.42 


88.39 

.. 

90.52 
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Table 2a - Energy Factor, J (m /sec ) 
Intermediate Flight Time Class 


Launch v ro 

Flight Time (days) 

IKSuKml 

550 

600 

650 

700 

800 

900 

0 

23.28 (20.3) 

12.75 (11.41) 

10.81 (10.5) 

10.41 (9.89) 

9.82 (9.26) 

9.71 (9.17) 

l 

17.76 

10.52 

10.20 


9.47 

9.24 

2 

13.34 

9.59 (9.04) 

9.62 

9.32 (8.89) 

9.13' 

8.87 

3 

9.96 

8.94 

9.05 

8.85 

8.84 

8.67 

4 

8.39 

8.37 (7.84) 

8.51 

8.41 (8.07) 

8.60 

8.61 

5 

7.63 

7.88 

8.00 

7.98 

8.39 


6 

7.13 

7.43 (7.00) 

7.50 

7.56 (7.29) 

8.20 


7 

6.70 

7.03 

7.03 

7.16 

8.00 


8 

6.53 

6.64 (6.25) 

6.59 

6.76 (6.35) 

7.81 



Table 2b - Propulsion Time Ratio, t /t 
Intermediate Flight Time Class ^ 





. Flight Time 

(days) 



550 

600 

650 

700 

800 

900 

0 

.6159 (1.0) 

.8650 (1.0) 

.8721 (1.0) 

.8748 (1.0) 

.8379 (1.0) 

.8165 (1.0) 

1 

.6913 

.8879 

.8617 


.7897 

.7949 

2 

.8003 

.8520 (1.0) 

.8513 

.8609 (1.0) 

.7440 

.7681 

3 

.9889 

.8357 

.8409 

.8512 

.7216 

.7492 

4 

.9014 

.8230 (1.0) 

.8305 

.8424 (1.0) 

.7105 

.7293 

5 

.8228 

.8108 

.8200 

.8338 

.7013 


6 

.7696 

.7977 (1.0) 

.8092 

.8251 (1.0) 

.6926 


7 

.7209 

.7824 

.7983 

.8164 

.6842 


8 

.7757 

.7428 (.843) 

. 7871 

.8075 (1.0) 

.6756 
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Table 2c - Mean Thrust Acceleration, a (10 m/sec ) 
Intermediate Flight Time Class 


3 

Flight Time (days) 

muSSSm! 

550 

600 

650 

700 

800 

900 

0 

8.918 

5.331 (4,70) 

4.698 

4.436 (4.05) 

4.118 (3.66) 

3.911 (3.44) 

1 

7.352 

4.781 

4.592 

n o 

4.165 

3.867 

2 

5.922 

4.661 (4.07) 

4.485 

4.231 (3.83) 

4.214 

3.853 

3 

4.604 

4.542 

4.379 

4.145 

4.209 

3.857 

4 

4.424 

4.430 (3.89) 

4.272 

4.062 (3.65) 

4.186 

3.897 

5 

4.418 

4.330 

4. 167 

3.978 

4.161 


6 

4.414 

4.240 (3.68) 

4.063 

3.893 (3.48) 

4.138 


7 

4.422 

4.162 

3.961 

3.807 

4.114 


8 

4.209 

4.153 (3.78) 

3.861 

3.720 (3.29) 

4.089 



28 , 




































Table 2d - Power to Mass Ratio, p /m (kw/kg) 
Intermediate Flight Time cfass° 


HB9 

Flight Time (days) 

550 

600 

650 

700 

800 

900 1 

0 

■SB 

m 

mm 

.0299 

mm 

mm 

1 

■Bi 

mBM 

■19 



mm 

2 

■9 

RIB 

■Hi 

.0287 

. 0287 


3 

mssm 

.0280 

.0283 

.0282 

.0284 

.0276 

4 

.0273 

.0273 

.0276 

.0276 

.0281 

.0275 

5 

.0263 

.0266 

.0269 

.0270 

.0278 


6 

.0256 

.0259 

.0262 

.0264 

.0276 


7 

.0251 

.0254 

.0255 

.0258 

.0273 


8 

.0248 

* 0248 

.0248 

.0252 

.0270 



Table 2e - Net Spacecraft Mass Ratio, m /m 
Intermediate Flight Time Class 


Launch v 

00 

Flight Time (days) 


550 

600 

650 

700 

800 

900 

0 

.5284 (.558) 

.6397 (.657) 

.6657 (.673) 

.6715 (.679) 

.6803 (.688) 

.6821 (.690) 

1 

.5812 

.6696 

.6743 


.6854 

.6892 

2 

.6317 

.6830 (.691) 

.6829 

.6876 (.694) 

.6905 

.6949 

3 

.6777 

.6931 

.6915 

. 6950 

.6950 

.6981 

4 

.7018 

.7020 (.710) 

.7000 

.7020 (.707) 

.6987 

.6989 

5 

.7142 

.7101 

.7085 

.7090 

.7021 


6 

.7229 

.7178 (.725) 

.7168 

.7160 (.720) 

.7052 


7 

.7305 

.7250 

.7251 

.7232 

.7084 


8 

.7339 

.7319 (.739) 

.7332 

. 7304 (. 734) 

.7117 



Table 2f - Optimum Jet Exhaust Speed, c (km/sec) 
Intermediate Flight Time Class 


Launch v 

CO 

Flight Time (days) 

[■mm 

550 

600 

650 

700 

800 

900 

0 

92.3 (114) 

115.7 (121) 

117.8 (125) 

122.9 (131.9) 

129.5 (143.5) 

134.9 (151.5) 

1 

98.5 

115.0 

117.4 


126.2 

133.3 

2 

106.6 

112.1 (122.8) 

117.0 

123.0 (131.5) 

123.0 

131.2 

3 

118.2 

111.1 

116.6 

122.7 

121.6 

129.5 

4 

111.0 

110.5 (117.5) 

116.1 

122.4 (128.6) 

120.9 

127.8 

5 

106.3 

109.9 

115.7 

122.1 

120.3 


6 

103.2 

109.3 (128.3) 

115.2 

121.7 (132.1) 

119.7 


7 

100.4 

108.5 

114.7 

121.3 

119.1 


8 

104.3 

106.1 (112.9) 

114.2 

121.0 (134.1) 

118.5 

— - 
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Table 3a - Energy Factor, J (m /sec ) 
Long Flight Time Class 


■inn 

■SB 

Flight Time (days) 

900 

1000 

1100 

1200 

1300 

0 

36.75 


(4.83) 

4.89 (4.68) 

4.29 

1 

29.09 



4.55 

3.70 

2 

22.56 

(4.013) 

4.40 

4.25 


3 . 

17.06 


4.22 

4.03 


4 

12.53 

2.34 (2.195) 

4.07 

3.84 


5 

8.88 

1.86 

3.92 

3.67 


6 

6.04 (5.78) 

1.56 (1.515) 

3.77 

3.51 


7 

3.94 

1.33 

3.63 

3.36 


8 

2.49 (2.455) 

1.13 (1.094) 

3.48 

3.20 



Table 3b - Propulsion Time Ratio, t /t 
Long Flight Time Class 15 


1 

flw93uSa338 

Flight Time (days) 

900 

1000 

1100 

1200 

1300 

0 • 

.5146 


(1.0) 

.9883 (.931) 

.9165 

1 

.5267 



.9488 

.8610 

2 

.5433 

(1.0) 

1.0 

.9129 


3 

.5640 


.9928 

, .8885 


4 

.5889 

1.0 (1.0) 

.9837 

.8772 


5 

.6180 

1.0 

.9718 

.8698 


6 

.6518 (.703) 

1.0 (1.0) 

.9610 

.8634 


7 

.6908 

1.0 

.9510 

.8574 


8 

.7350 (.702) 

.9804 (1.0) 

.9416 

.8514 
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Table 3c - Mean Thrust Acceleration, a (10 m/sec ) 
Long Flight Time Class 



Flight Time (days) 

900 

1000 

1100 

1200 

1300 

0 

9.58 


(2.26) 

2.19 (2.20) 

2.04 

1 

8.43 



2.15 

1.96 

2 

7.31 

(2.15) 

2.15 

2.12 


3 

6.24 


2.11 

2.09 


4 

5.23 

1.64 (1.59) 

2.09 

2.06 


5 

4.30 

1.47 

2.06 

2.02 


6 

3.45 (3.25) 

1.34 (1.32) 

2.03 

1.98 


7 

2.71 

1.24 

2.00 

1.94 


8 

2.09 (2.12) 

1.16 (1.125) 

1.97 

1.91 
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Table 3d - Power to Mass Ratio, p /m (kw/kg) 
Long Flight Time Class 


Launch v 

00 

(km/sec) 

Flight Time (days) 

900 

1000 

1100 

1200 

1300 S 

0 

msm 

mu 

mam 

mm 

mm 

1 

H9 


mMm 

WBm 


2 

KB 



WSSM- 

BB 

3 

wBm 


BB 

.0202 


4 

.0338 


B 

.0198 


5 

.0292 


RSI 

.0194 


6 

. 0247 

.0130 

■■ 

.0190 


7 

.0203 

.0121 

SBB 

.0186 


8 

.0164 

.0112 

.0189 

.0182 

■■ 


Table 3e - Net Spacecraft Mass Ratio, m /m 
Long Flight Time Class ; 


Launch v 

09 

(km/sec) 

Flight Time (days) 

900 

1000 

1100 

1200 

1300 

i 

0 

.4313 

j 

(.770) 

.7690 (.773) 

mm i 

1 

.4837 



.7768 


2 

.5366 

(.789) 

.7802 

.7836 

■ 

3 

.5896 


.7844 

.7891 


4 

.6423 

.8372 (.841) 

.7882 

.7937 


5 

.6941 

.8540 

.7918 

.7981 


6 

.7440 (.749) 

.8658 (.867) 

.7955 

.8023 


7 

.7906 

.8760 

.7993 

.8065 

B B 

8 

.8318 (.831) 

.8851 (.887) 

.8032 

.8107 



Table 3f - Optimum Jet Exhaust Speed, c (km/sec) 
Long Flight Time Class 


Launch v 

00 

(km/sec) 

Flight Time (days) 

900 

1000 

1100 

1200 

1300 

0 

105.4 


(169.5) 

179.0 (174.3) 

180.1 

1 

108.0 



175.9 

174.6 

2 

111.1 

(163.6) 

173.1 

172.7 


3 

114.5 


172.6 

170.5 


4 

118.3 

167.9 (157.5) 

171.9 

169.5 


5 

122.4 

167.4 

171.0 

168.9 


6 

126.8 (131.0) 

167.6 (164.9) 

170.2 

168.4 


7 

131.6 

167.8 

169.5 

167.9 


8 

136.8 (118.4) 

166.5 (167.5) 

168.9 

167.5 
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Table 4a - Energy Factor, J (m /sec ) 
Extra Long Flight Time Class 


Flight Time (days 


1300 

1350 

1400 

25.13 

19.44 



14.58 

10.51 

7.20 

4.59 


1.62 

2.63 


1.30 (1.276) 

1.29 

.41 

1.10 

.51 

.29 

.93 



Table 4b - Propulsion Time Ratio, t /t^ 
Extra Long Flight Time Class P 




--4.2 

Table 4c - Mean Thrust Acceleration, a (10 m/sec ) 
Extra Long Flight Time Class 


Flight Time (days 



7.28 



6.42 



5.51 



4.61 



3.73 



2.88 


1.16 

2.09 


1.04 (1.027) 

1.39 

.59 

.95 

.82 

.50 

• 

00 

00 
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Table 4d - Power to Mass Ratio, p /m (kw/kg) 
Extra Long Flight Time cfass° 


Launch 


Flight Time (days) 

(km/sec) 

1300 

| 1350 

1400 

0 

.0448 

mmsmm. 


1 




2 

.0362 



3 

.0316 



4 

.0268 



5 

.0219 


.0134 

6 

.0170 


.0120' 

7 

.0121 


.0111 

8 

.0077 


.0102 

r 

Table 4e - Net Spacecraft Mass Ratio, m /m 


Extra Long Flight Time Class 

u u 

Launch v 

00 

Flight Time (days) 


1300 

1350 

1400 

0 

.5151 



l 

.5658 



2 

.6175 



3 

.6698 



4 

.7224 



5 

.7751 


.8636 

6 

.8273 


.8774 (.878) 

7 

.8777 

.9304 

.8871 

8 

.9225 

.9413 

.8957 


Table 4f - Optimum Jet Exhaust Speed, c (km/sec) 
Extra Long Flight Time Class 



Flight Time (days) 

1300 

1350 

1400 

0 




1 

Jiiil- *y,-. 



2 




3 

^ H 



4 

128.5 



5 

133.9 


200.4 

6 

140.8 


200.1 (195) 

7 

149.7 

198.5 

200.3 

8 

159.3 

198.5 

200.6 


























Table 5 - Performance Parameters for 700 Day Encke Rendezvous Mission 

Multi-parameter Independent Mode 
Arrival on October 17, 1980 


Launch v 

OO 

km/ sec 

J 

2 , 3 

m /sec 

t /t. 
P f 

a 

-4 2 

10 m/sec 

p /m 
o o 

kw/kg 

V 

n 

c 

km/sec 

m /p 
n o 

kg/kw 

0 

— 

10.355 

1.0 

4.1379 

■ — 1 " 

.0059075 

.37285 

27.946 

63.115 

1 

9.890 

1.0 

4.0438 

.0057949 

.37938 

27.776 

65.467 

2 

9.441 

1.0 

3.9510 

.0056625 

.38462 

27.493 

67.924 

3 

9.000 

1.0 

3.8576 

.0055248 

.38970 

27.179 

70.536 

4 

8.564 

1.0 

3.7630 

.0053844 

.39485 

26.849 

73.332 

5 

8.134 

1.0 

3.6673 

.0052425 

.40014 

26.508 

76.327 

6 

7.710 

1.0 

3.5704 

.0050992 

.40562 

26.159 

79.545 

7 

7.292 

1.0 

3.4723 

.0049551 

.41130 

25.802 

83.006 

8 

6.881 

1.0 

3.3731 

.0048101 

.41721 

25.438 

86.736 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 


100 


Energy 

Factor, 

J 

(m^/sec^) 



400 


600 


800 


1000 


1200 


1400 


Flight Time (days) 

Figure 2 - Variation of energy factor with flight time. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 



400 600 800 1000 1200 1400 

Flight Time (days) 

Figure 3 - Variation of mean thrust acceleration with flight time. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 



r nttn n ttrt lllil l Hi t mt t t III 1 1 1111)1 ttttttttt tt 1 I 111 1 1.1. HU I 1 1 

400 600 800 1000 1200 1400 

Flight Time (days) 

Figure 4 - Variation of propulsion time ratio with flight time. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 



* 3330 3530 3730 3930 4130 4330 4530 

Julian Date (-2440000) 

Figure 5 - NEP engine switch function time histories. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 


Energy- 

Factor, 

J 

(m 2 /sec^) 



0 2 4 6 8 

Launch (km/ sec) 


Figure 6- Variation of energy factor with launch excess speed. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 


Mean 

Thrust 

Acceleration, 

a 

(m/sec 2 ) 



Figure 7 - Variation of mean thrust acceleration with launch excess speed. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 


Propulsion 

Time 

Ratio, 




Launch v (km/sec) 

00 

Figure 8 - Variation of propulsion time ratio with launch excess speed. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrive on October 17, 1980 
Intermediate Flight Time Class 


Energy 

Factor, 

J 

(m 2 /sec 3 ) 



0 10 20 30 40 


Specific Propulsion System Mass, & (kg/kw) 


Figure 9 - Variation of energy factor with specific propulsion system mass. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 
Intermediate Flight Time Class 


Mean 

Thrust 

Acceleration 

a 

(m/sec ) 



0 10 20 30 40 


Specific Propulsion System Mass, a (kg/kw) 


Figure 10 - Variation of mean thrust acceleration with specific propulsion 
system mass. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
Arrival on October 17, 1980 
Intermediate Flight Time Class 


Propulsion 

Time 

Ratio, 




Figure 11 


Variation of propulsion time ratio with specific propulsion system 
mass. 
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1980 NEP ENCKE RENDEZVOUS MISSION 


Specified t , t_, m and v 's 
o f o 00 

Optimized c 



Figure 12 - Typical behavior of net spacecraft mass as a function of 
power level. 
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1980 NEP ENCKE RENDEZVOUS MISSION 
700 Day Flight Time Arriving on October 17, 1980 
Titan III D/Centaur Launch Vehicle 
a. = 5 kg/kw 



Figure 13 - Performance Parameters for Encke Mission at Optimum Power 
Level. 
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Figure 14 


1980 NEP ENCKE RENDEZVOUS MISSION 
700 Day Flight Time Arriving on October 17, 1980 
Tital III D/Centaur Launch Vehicle 
a = 5 kg/kw 



0 2 4 6 8 


Launch Excess Speed (km/sec) 

- Performance Parameters for Encke Mission at Constrained Power 
Levels using Multi-Parameter Independent Mode. 
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